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The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, 
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unresolved ring. General relativity predicts that embedded within this image lies a thin “photon ring,” which is 
composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around 
the black hole. The subrings approach the edge of the black hole “shadow,” becoming exponentially narrower but 
weaker with increasing orbit number, with seemingly negligible contributions from high-order subrings. Here, we 
show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures 
offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only 


a sparse interferometric array. 


INTRODUCTION 

The Event Horizon Telescope (EHT) Collaboration has recently 
published images of the supermassive black hole in M87 using very 
long baseline interferometry (VLBI) at 1.3-mm wavelength (1-6). 
These images reveal a bright ring of emission with a diameter of 
approximately 40 uas. However, while the diameter of this ring is 
resolved by the EHT, its thickness and detailed substructure are not. 
Here, we show that general relativity predicts an intricate substructure 
in the ring that presents distinctive signatures for interferometric 
measurements. These signatures offer a promising approach for 
precisely determining the mass and spin of black holes and for testing 
general relativity using sparse interferometers, such as an extension 
of the EHT to space. 

Neglecting opacity, a telescope with perfect resolution directed 
at a black hole observes an infinite number of nested images of the 
universe. These images arise from photons that differ by the number 
n of half-orbits that they complete around the black hole on the way 
from their source to the detector. Each such image is thus an in- 
creasingly delayed and demagnified snapshot of the universe as seen 
from the black hole. In an astrophysical setting, this self-similar 
sequence of relativistic images is dominated by the luminous matter 
surrounding the black hole and produces in its image a feature 
known as the “photon ring” of the black hole (7-10). The leading 
(n = 1) subring appears as a sharp, bright feature in ray-traced images 
from many general relativistic magnetohydrodynamic (GRMHD) 
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simulations (see Fig. 1). Successive subrings have exponentially sharper 
profiles and asymptotically approach the boundary of the black hole 
“shadow.” For large n, these profiles mirror the leading subring in a 
manner that universally depends on the spacetime geometry, with 
the ratio of successive subring flux densities determined by Lyapunov 
exponents that characterize the instability of bound photon orbits. 
Hence, measuring the size, shape, and thickness of the subrings 
would provide new and powerful probes of a black hole spacetime. 

Both GRMHD simulations and analytic estimates suggest that 
the photon ring should provide only ~10% of the total image flux 
density. This dimness may appear to preclude observations of the 
photon ring and its substructure, which is dimmer still. However, 
interferometric measurements are sensitive to more than just overall 
flux: They also natively filter images by their spatial wave numbers 
and therefore naturally isolate contributions from individual photon 
subrings. Sufficiently long baselines also resolve out diffuse flux in 
an image and are thus dominated by power from the photon ring. 
Hence, although sharp elements of the photon ring produce a negli- 
gible contribution to the total flux in an image, they can still provide 
a pronounced, dominant signal on long baselines. 

In this Research Article, we explore the photon ring’s theoretical 
underpinnings and show that, unexpectedly, precise measurements 
of the photon ring and even its subrings are feasible using inter- 
ferometry. We describe the shell of bound photon orbits of a Kerr 
black hole and its relation to the photon ring. We also present a 
decomposition of the photon ring into subrings indexed by half-orbit 
number and derive their self-similar structure, which is universally 
governed by Lyapunov exponents that characterize orbital instability. 
We then derive generically expected interferometric signatures of 
the photon ring. We show that its subrings produce a cascade of 
damped oscillations on progressively longer baselines, with the 
visibility of each subring conveying precise information about its 
diameter, width, and angular profile. Last, we discuss observational 
prospects for detecting these signatures with extensions of the EHT. 
In particular, we highlight the possibility of detecting the leading 
n = 1 subring using a station in low Earth orbit, the = 2 subring 
using a station on the Moon, and the n = 3 subring using a station 
in the Sun-Earth L, orbit. 
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Fig. 1. Time-averaged image of a GRMHD simulation of M87. This model has param- 
eters chosen to be consistent with the 2017 EHT data and corresponds to the high 
magnetic flux “magnetically arrested disk” accretion state. It has M = 6.2 x 10°Mo, 
a/M=0.94, Bobs = 163°, Fhigh = 10, and mass accretion rate matching the 1.3-mm flux 
density (5). The spin axis points left when projected onto the image. The time average 
was performed over 100 snapshots produced from uniformly spaced GRMHD fluid 
samples over a time range of 1000M (approximately 1 year). Although visually 
prominent, the thin, bright ring contains only ~20% of the total image flux density. 


RESULTS 

Photon shell and photon ring 

This section describes the shell of unstable bound photon orbits 
surrounding a black hole, its lensed photon ring image, the photon 
subrings labeled by half-orbit number, and the angle-dependent 
Lyapunov exponents that govern the subring brightness ratio asymmetry. 
Previous treatments of these structures include (8) and (10-13). 
Photon shell 

The photon shell, illustrated in Fig. 2, is the region of a black hole 
spacetime containing bound null geodesics or “bound orbits” that 
neither escape to infinity nor fall across the event horizon. For 
Schwarzschild, the photon shell is the two-dimensional sphere at 
r = 3M and any 8, 9, and t. For Kerr, this two-dimensional sphere 
fattens to a three-dimensional spherical shell. It is best described 
using Boyer-Lindquist coordinates, in which the metric of a Kerr 
black hole of mass M and angular momentum J = aM (with 0 < 
a<M)is 


ds? = -2 (at - asin? @do) + Žar + Edo? 


+2 
+ nO? + a”) do — adt]? 


(1A) 
A = r?-2Mrt+a’, E = r*+a’cos*0 (1B) 
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Fig. 2. Photon shell and photon ring of a Kerr black hole. (A) Cross section of the 
photon shell in the (r, 6) plane in Boyer-Lindquist coordinates. The black hole spin 
is a/M = 0.94, directed vertically, and the color varies with r. The intersection of an 
observer's line of sight with the photon shell boundaries at r = r} determines the 
visible subregion of the photon shell. (B to D) Photon ring on the screen of an observer 
at varying inclinations Oop; relative to the spin axis, whose projection onto the plane 
perpendicular to the line of sight is depicted by the (left-pointing) arrow. The center 
of the photon ring has a displacement from the origin that increases with spin. 
The color coding on the ring denotes the matching radius on the shell from which the 
photon emanated. The photon shell r < r < r* is only visible in its entirety to 
the edge-on (Oops = 90°) observer. The face-on (Obs = 0°) observer only receives 
photons from the white r = nr orbit. The Oop; = 17° observer sees the portion of 
the shell delineated by the dashed lines. 


These coordinates have the special property that each bound or- 
bit lies at some fixed value of r in the range 


(2A) 


ri = 2M|1 + cos (Zarccos (+ +))| (2B) 


Every point in the equatorial annulus r’ < r < r‘,@=7n/2,has 
a unique bound orbit passing through it. On the boundaries r = r}, 
the orbits reside entirely in the equatorial plane. At generic points, 
on the other hand, they oscillate in the 6 direction between polar angles 


0. = arccos (+u,) (3) 
where 
Uy = impl- +3M°r-2a°M 
+ MA (2r? -3Mr” + a? M) | (4) 


We will refer to one such complete oscillation (e. g., from 80- back to 
itself) as one orbit, since the photon typically returns to a point near, but 
not identical to (since the azimuthal angle ọ also shifts), its initial position. 
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To summarize, the photon shell is the spacetime region 


ro<r<ri, 0- <0<0, 0<0<2n (5) 
(depicted in Fig. 2) for all times — œ < t < œ. 
The bound orbit at radius r has the energy-rescaled angular 


momentum 


fe M(r?—a?)-rA 
~ alr- M) 


The inner circular equatorial orbit at r” is prograde, while the outer 
one at r” is retrograde: €(r) 2 0. The overall direction of the orbits 
reverses at the intermediate value rX for which / vanishes. At that 
radius, [0-, 0+] equals [0, 2], and the orbits can pass over the poles. 

The bound geodesics are unstable in the sense that, if perturbed 
slightly, they either fall into the black hole or escape to infinity where 
they can reach a telescope. The observed photon ring image arises 
from photons traveling on such “nearly bound” geodesics. Consider 
two geodesics, one of which is bound, with the other initially differing 
only by an infinitesimal radial separation ôro. The equation of geodesic 
deviation shows that, after n half-orbits between 04, their separation 
grows to 


(6) 


õrn = e" ôro 


(7) 


Here, the so-called Lyapunov exponent y is a function on the space 
of bound orbits given by (see the Supplementary Materials) 


_A r2 MrA 
(r- M) 


(8) 


1 dt 
Jaa - #7) (ust? -u 


A closely related formula appears in (14). Hence, the nearly bound 
geodesic will typically cross the equatorial plane a number of times 
of order 


Ory 
oro 


until ôr, >> õro, when the geodesic is well separated from the bound 
orbit and it shoots off to infinity (or crosses the event horizon if ôro < 0). 
These Lyapunov exponents are central and potentially observable 
quantities that characterize the geometry of the Kerr photon shell. 
Photon ring and subrings 

The photon ring is the image on the observer screen produced by 
photons on nearly bound geodesics (7). In the limit in which the 
photons become fully bound, it may be shown that their images ap- 
proach a closed curve C, given by 


ne qin (9) 


p= D7! \a2(cos”O obs — u,u_) +27 (10A) 


= _ £ 
Pp = arccos( Den, ET) (10B) 


where (p, Pp) are dimensionless polar coordinates on the observer 
screen, while (D, obs) denote the observer’s distance and inclination 
from the Kerr spin axis, respectively. We can view C, as parameterized 
by the shell radius r” < r < r! from which the photon originated. 
For each value of r, Eq. 10B has two solutions for @p in the range 0 
< p < 27, so each radius in the photon shell appears at two posi- 
tions on Cy. A notable consequence of Eqs. 10A and 10B is that for 
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Oobs # 0, both 7 and p, and hence @p, are functions only of r, Bobs, 
and D. Hence, a measurement at a specific angle pp along the ring 
probes a specific radius r of the Kerr geometry and not, as might 
have been expected, a specific angle around the black hole! 

Astrophysically observed photon intensities I,ing(p, Pp) at the screen 
can be computed by backward ray tracing. One follows the null geo- 
desics from the observer screen back into the Kerr spacetime, 
integrating the Doppler-shifted strength J of matter sources along 
the geodesic, with attenuation factors accounting for the optical 
depth. Scattering effects are negligible because the expected plasma 
frequency and electron gyroradius are in the megahertz range, several 
orders of magnitude below the observing frequencies that we con- 
sider. For the images in this paper, we used ipole (15). A light ray 
aimed exactly at the curve C, is captured by the photon shell and 
(unstably) orbits the black hole forever. Those aimed inside C, fall 
into the black hole, while those aimed outside escape to infinity. 
Therefore, Cy is the edge of the black hole shadow. 

If we shoot a light ray very near, a distance dp from the shadow 
edge at pa it will circle many times through the emission region before 
falling into the black hole or escaping to infinity. The affine length 
of the ray and its number of half-orbits accordingly diverge as 5p — 0 


(11) 


This follows from Eq. 9 together with a computed relation be- 
tween dp and dro. For optically thin matter distributions, Eq. 11 implies 
a mild divergence in the observed ring intensity ling ~ n as the 
shadow edge is approached, since a light ray that completes n half- 
orbits through the emission region can collect ~n times more photons 
along its path. The photon ring is then the bump in the photon in- 
tensity containing this logarithmic divergence at the shadow edge. 
Although the divergence is cut off by a finite optical depth, this notable 
feature remains visually prominent in many ray-traced images of 
GRMHD simulations, as in Fig. 1. 

The photon ring can be subdivided into subrings arising from 
photons that have completed n half-orbits between their source and 
the screen. This definition for the photon ring agrees with that in (16) 
but differs from the later usage in (9) and (10) by the inclusion of 
the n = 1 and 2 contributions. These low n contributions fully account 
for the thin ring image visible in Fig. 1. To orbit at least n/2 times 
around the black hole, the photon must be aimed within an expo- 
nentially narrowing window 


5p n 

Pc 
around the shadow edge. Hence, the subrings occupy a sequence of 
exponentially nested intervals centered around C}. 

Each subring consists of photons lensed toward the observer screen 
after having been collected by the photon shell from anywhere in 
the universe. Hence, in an idealized setting with no absorption, each 
subring contains a separate, exponentially demagnified image of the 
entire universe, with each subsequent subring capturing the visible 
universe at an earlier time. Together, the set of subrings are akin to 
the frames of a movie, capturing the history of the visible universe 
as seen from the black hole. In an astrophysical setting, these images 
are dominated by the luminous matter around the black hole. For a 
black hole surrounded by a uniform distribution extending over the 
poles, the contributions made by each subring to the total intensity 
profile cannot be told apart, and the individual subrings cannot be 


-yn 


xe (12) 
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distinguished on the image. However, for a realistic disk or jet with 
emission peaked in a conical region, the subrings are visibly distinct: 
the nth subring is approximately a smooth peak of width e ””. Summing 
these smooth peaks, like layers in a tiered wedding cake (see Fig. 3), 
reproduces the leading logarithmic divergence in the intensity (Eq. 11). 

The photons comprising successive subrings for the same angle 
@p traverse essentially the same orbits and hence encounter the same 
matter distribution around the black hole. Apart from source varia- 
tions on the time scale of an orbit, intensities of the nth and (n + 1)th 
subring differ only because they correspond to windows whose 
widths 5p, and 5p, + ı differ by a factor of e”. Hence, for large enough 
n, the intensities are related by 


Iing Pe + 5P, Pp) S Iing(Pe + e” òp, Pp) (13) 
We therefore find the angle-dependent subring flux ratio 
Fe 
mom eT (14) 


Equations 13 and 14 are matter-independent predictions for the 
photon ring structure that involve only general relativity. The pre- 
diction holds only for “large enough” n: At small n, there are non- 
universal matter-dependent effects from photons that do not traverse 
exactly the same region around the black hole. Insight into when n 
is large enough might be obtained from GRMHD simulations. 

Since the exponent y depends on a, Bobs, and Pp (see the Supple- 
mentary Materials), the flux ratio asymmetry in Eq. 14 provides a 
new method for determination of the spin. For Schwarzschild, y = n 
(8), corresponding to a demagnification factor of e* + 4%. For a 
black hole of maximal spin a/M = 1 viewed from an inclination Oops = 
17° [as estimated for M87; (17)], the factor e™” is as large as 13% on 
the part of the ring where the black hole spins toward the observer. 
Although Eq. 12 breaks down for n = 0, this suppression factor 
suggests that the leading n = 1 subring should provide ~10% of the 
total luminosity, in order-of-magnitude agreement with GRMHD 
simulations. 


Interferometric signatures of a photon ring 

This section explores the response of an interferometer to the 
photon ring through a series of increasingly refined geometrical 
models. We first review the interferometric signatures of an infini- 
tesimally thin, uniform, and circular ring. We then extend this 
treatment to include rings with nonuniform brightness, nonzero 
thickness, and noncircular structure. We conclude this section 
by discussing specific features expected for the photon ring and 
its subrings. 

Visibilities for a thin, uniform, and circular ring 

Each baseline joining two elements of an interferometer samples a 
complex visibility V(u), which corresponds to a single Fourier com- 
ponent of the sky image I(x) (18) 


Vu) = fi(x) e7™* d?x (15) 


Here, u is the dimensionless vector baseline projected orthogonal 
to the line of sight and measured in units of the observation wave- 
length A, while x is a dimensionless image coordinate measured 
in radians. 
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Fig. 3. Image cross sections of a photon ring and its subrings. (A) Brightness 
cross sections for the time-averaged GRMHD image shown in Fig. 1. The blue/red 
curves show cross sections perpendicular/parallel to the projected spin axis. (B and 
C) Decomposition of the left perpendicular peak and the right parallel peak into 
subrings indexed by the number n of photon half-orbits executed between turn- 
ing points (Eq. 3) in the polar motion. Similar results are also seen in image cross 
sections of simple geometrical models (10). 


In terms of polar coordinates (p, Pp) on the observer screen (Eqs. 
10A and 10B), the image and corresponding visibility function of an 
infinitesimally thin, uniform, and circular ring are 


d 
z) 
V(u, Pu) = Jo(ndu) 


where d is the ring diameter in radians and the image is normalized 
to have a total flux density of unity, V(0) = 1. Jm denotes the mth Bessel 
function of the first kind, which admits the asymptotic expansion 


In(adu) = 12 cos [= (au — 2m+1)] 


valid for ndu > m?. Hence, V(u) is a weakly damped pure frequency 
with period Au = 2/d inside an envelope that falls as 1/ vu. 


I(p,@p) = 13(0 (16A) 


(16B) 


(17) 
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Fig. 4. Universal interferometric signatures of a photon ring. (A to D) Visibility amplitudes of (A and B) the time-averaged GRMHD simulation shown in Fig. 1 and 
(C and D) a GRMHD snapshot (see the Supplementary Materials). Amplitudes are shown for baselines perpendicular (red) and parallel (blue) to the black hole spin axis. 
While short baselines (left of the vertical dotted lines) display complex structure reflecting astrophysical features of the image such as emission from the disk and jet, 


longer baselines are dominated by the universal interferometric signatures of the photon ring. A simple model | V(u)| = | a+ cos (ndu) + a- sin (zdu) | (du) Ves is 
overplotted (black dashed curves), with parameters determined independently along the two axes. The periodicities encode the ring diameters along each axis and 
hence M/D for the black hole; their difference provides an estimate of the black hole spin and inclination. The parameters a carry information about the angular bright- 


ness distribution (and hence spin and inclination). The dashed green curve yo? 


shows the expected envelope for an infinitesimally thin ring, while the solid green curve 


ue". shows the fitted envelope that carries information about the ring thickness. On even longer baselines (B and D), the dominant visibility contributions arise from 
subrings with increasingly higher n. The universal features are more prominent in the time-averaged image, whose ring is dominated by smaller mode numbers m and 


which has less small-scale power outside the photon ring. 


Visibilities for a nonuniform ring 
The image of a thin ring with nonuniform brightness in p decom- 
poses into a sum over angular Fourier modes 


psp) = F(p -4) p> Bme!” 


where B-m = Bh since the image is real. The total image flux density 
is given by fo > 0. 
The corresponding visibility function is 


(18) 


co 


Vee Pu) = YY BmJmn(mdu) e0 (19) 
Using Eq. 17, for long baselines, we may approximate 
04(@,) cos (tdu) +a_(@,,) sin (ndu) 
V(u, Qu) S 
(u, Pu) RET 
alat È pet bora] 
i (20) 


Thus, for sufficiently long baselines, the radial visibility function 
of a nonuniform thin ring is determined by a single pair of weakly 
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damped, orthogonal modes a.(@,,). Their envelope still fall as 
| V(u)| ~ 1/vu, and the modes have a common period of Au = 2/d 
in complex visibilities (or Au = 1/d in visibility amplitudes). The 
angular spectrum of the image {Bm} is easily retrieved from the 
angular spectrum of the visibilities (see the Supplementary Materials). 
Visibilities for a thick ring 

Baselines of length u 2 1/L are required to resolve image features of 
size SL. Hence, the visibility function of any ring with diameter d 
and thickness w « d has two asymptotic regimes 


M:i kuci M4 «j< uv (21) 

Baselines in regime (I) resolve the diameter of the ring but not its 
thickness, while longer baselines in regime (II) resolve both. Hence, 
the visibility function in regime (I) behaves like that of a thin ring 
(a damped periodicity with envelope | V(u)| ~ 1/Vu), while the 
envelope of the visibility function in regime (II) is sensitive to the 
radial profile of the ring. In general, the visibility of any smooth ring 
decays exponentially in regime (II), although images with discontinuous 
derivatives, such as a uniform disk or annulus, can have slower, 
power-law falloffs. 
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Fig. 5. Prospects for observing a photon ring. (A) Schematic showing visibility 
amplitude as a function of baseline length for a photon ring with d = 40 pas and 
flux density comparable to M87. The black, cyan, and magenta visibilities correspond 
to photons with half-orbit numbers n= 1, 2, and 3. (B) Frequency-dependent range 
of Earth baselines and representative Earth-space baselines. Earth-space baselines 
shown are the longest baselines for an orbiter in low Earth orbit (LEO), in medium 
Earth orbit (MEO), in geostationary orbit (GEO), on the Moon, and at the second 


Sun-Earth Lagrange point (L2). 


The validity of the approximation of Eq. 17 in regime (I) depends 
on the amount of power at high values of m. Specifically, it requires 
Mmax < Vnd/w. Under this condition, | V(u) | has <d/w periods in 
regime (I). 

Visibilities for a noncircular ring 

Although the photon ring is nearly circular for all black hole spins 
and inclinations, the primary interferometric signatures discussed 
thus far do not require an image with perfectly circular structure. 
For instance, if an image is stretched, I(x, y) > I(x, y) = I(ax, by); 
then, its visibility function is correspondingly compressed, V(u, v) > 
V' (u,v) = | ab| ~'V(u/a, v/b). Thus, the visibility profiles of a stretched 
ring share the properties and asymptotic expansions derived for a 
circular ring (e. g., Eq. 20), except that the radial periodicities become a 
function of position angle. To leading order in the asymmetry 1 — a/b, 
the diameter corresponding to a damped radial periodicity in the 
visibility domain matches that of the stretched ring along the base- 
line’s position angle. For the black hole in M87, the asymmetry is 
expected to be a few percent at most, even for a maximally rotating 
black hole (see fig. S2). 

Visibilities of the photon subrings 

As discussed earlier, the photon ring decomposes into subrings 
labeled by the photon half-orbit number n. According to Eq. 14, the 
width of the radial intensity profile produced by the nth subring is 
Wn ~ woe '”, while the brightness remains approximately constant 
with n (until some nmax determined by the optical depth). Each sub- 
ring thus contributes a periodically modulated visibility, V,(u) ~ 
w,/\u, which falls more steeply for baselines u > 1/w,. Hence, the 
nth subring dominates the signal in the regime 


1 1 
Wnai Ku x Wa 


(22) 


This implies that the totality of subring contributions has an 
envelope defined by this turnover behavior 
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re w 1 
V(u)~ È E~ (23) 


Wr<l/u á 

Together, the subrings then form a cascade of damped oscillations 
on progressively longer baselines, each dominated by the image of a 
single subring and conveying precise information about its diameter, 
thickness, and angular profile. Figure 4 displays visibilities of the 
time-averaged GRMHD image in Fig. 1, which exhibit the expected 
damped periodicity, as well as clear contributions on long baselines from 
distinct subrings (see Fig. 5 for a schematic illustration of this cascade). 


DISCUSSION 

Detection of the photon ring’s universal interferometric signatures 
requires measurements on longer baselines, with finer angular reso- 
lution than those currently available to the EHT. This extension can 
be achieved either by observing at higher frequencies or on longer 
physical baselines via space VLBI. For reference, ALMA (Atacama 
Large Millimeter/submillimeter Array) currently observes up to 
950 GHz (ALMA band 10) with higher frequencies (up to 1.53 THz) 
possible in the future (19, 20). Figure 5 shows baseline lengths for a 
variety of array configurations and observing frequencies. 

For the EHT to observe the photon ring, it must also achieve 
sufficient sensitivity to detect its visibilities. For both Sgr A* and M87, 
the horizon-scale emission has a total flux density of Fiot ~ 1 Jy at 
A~1 mm (21, 22), with a fraction fring ~ 10% expected to come from 
the photon ring. The expected amplitude of the photon ring on long 
interferometric baselines is thus 


-3/2 aija fF 
a en ay la o (a a (5) (a) (24) 


For comparison, a baseline from ALMA to a 4-m orbiter with 
32 GHz of averaged bandwidth and a 10-min coherent integration 
would have a thermal noise of 6959 = 3 mJy at 950 GHz and O69 = 
1.3 mJy at 690 GHz. For baselines from ALMA to a 10-m orbiter, 
such as the proposed Millimetron mission for L3 (23, 24), the thermal 
noise would be 959 ~ 1 mJy and o¢99 ~ 0.5 mJy. Another possibility 
would be to place a VLBI station on or orbiting the Moon, which 
could sample many periods of the n = 2 regime of M87 at current 
EHT observing frequencies. A 10-m dish on the Moon could 
achieve o ~ 0.1 mJy on baselines to ALMA with coherent integra- 
tions of 10 min and a bandwidth of 32 GHz. 

These sensitivities could be substantially improved via simultaneous 
multifrequency observations. In addition to having more sensitive 
receivers and longer coherence times, lower frequencies give corre- 
spondingly shorter baselines and thus increased interferometric power 
from the photon ring. Phase calibration with lower frequencies could then 
allow substantially longer integration times at higher frequencies (see 
the Supplementary Materials for additional details and discussion). 

Interferometric signatures of the photon ring are most prominent 
when the image has little small-scale power outside the ring and 
when the ring has a smooth angular profile dominated by low mode 
numbers m. Both of these conditions are met in time-averaged images 
of black hole accretion flows, such as in Fig. 1. Because visibilities of 
a time-averaged image are equal to time-averaged visibilities of a 
variable image, developing capabilities for long, coherent averaging 
could significantly improve the prospects for unambiguous detec- 
tion and characterization of the photon ring. 
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In summary, precise measurements of the size, shape, thickness, and 
angular profile of the nth photon subring of M87 and Sgr A* may be 
feasible for n = 1 using a high-frequency ground array or low Earth orbits, 
for n = 2 with a station on the Moon and for n = 3 witha station in L}. 


SUPPLEMENTARY MATERIALS 

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/ 
content/full/6/12/eaaz1310/DC1 

Supplementary Materials and Methods 

Fig. S1. Lyapunov exponent as a function of image angle. 

Fig. S2. Black hole shadow diameter and asymmetry as a function of spin and inclination. 
Fig. S3. Snapshot images and visibilities. 
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